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ABSTRACT 

The supernova (SN) PTFlliqb was initially classified as a Type Iln event caught 
very early after explosion. It showed narrow Wolf-Rayet (WR) spectral features on 
day 2 (as in SN 1998S and SN 2013cu), but the narrow emission weakened quickly 
and the spectrum morphed to resemble those of Types II-L and II-P. At late times, Ha 
emission exhibited a complex, multipeaked profile reminiscent of SN 1998S. In terms of 
spectroscopic evolution, we find that PTFlliqb was a near twin of SN 1998S, although 
with a factor of ^ 2-4 weaker interaction with circumstellar material (CSM) at early 
times, and stronger CSM interaction at late times. We interpret the spectral changes as 
caused by early interaction with asymmetric CSM that is quickly (by day 20) enveloped 
by the expanding SN ejecta photosphere, but then revealed again after the end of the 
plateau when the photosphere recedes. The light curve can be matched with a simple 
model for weak CSM interaction (with a mass-loss rate of roughly 10“^ Mq yr“^) 
added to the light curve of a normal SN II-P; the relatively weak CSM interaction 
allows the plateau to be seen more clearly than in other SNe Iln. This plateau requires 
that the progenitor had an extended hydrogen envelope like a red supergiant at the 
moment that it exploded, consistent with the slow progenitor wind speed (< 80 km 
s“^) indicated by narrow Ha emission. The cool supergiant progenitor is significant 
because PTFlliqb showed WR features in its early spectrum — meaning that the 
presence of such WR features in an early SN spectrum does not necessarily indicate a 
WR-like progenitor. While the late-time Ha profile was multipeaked and asymmetric 
like that of SN 1998S, PTFlliqb’s Ha developed a stronger redshifted peak, so in 
this case the asymmetry cannot be blamed on dust obscuration. Instead, azimuthal 
asymmetry due to binary interaction is likely. Overall, PTFlliqb bridges SNe Hn with 
weaker pre-SN mass loss seen in SNe H-L and H-P, implying a continuum between 
these types. It hints at episodic pre-SN mass loss on a wider scale than previously 
recognised. 

Key words: circumstellar matter — stars: evolution — stars: winds, outflows — 
supernovae: general — supernovae: individual (PTFlliqb) 
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1 INTRODUCTION 

Recent years have seen the growing recognition that episodic 
mass loss may be a critical ingredient in the evolution of 
massive stars. Perhaps the most vivid demonstration is the 
class of Ty pe Iln supernovae (SNe), which h ave narrow emis¬ 
sion lines ((Schlegel]ll990l : lFiIiDpentolll997l ~l arising from ex¬ 
tremely dense circumstellar material (CSM) that is close to 
the star and must have been ejected in o nly the f ew ye ars or 
decades preceding the final explosion; see lSmi t3 (l2014h for a 
general review of massive-star mass loss and its connection 
to SN diversity. 

The synchronisation of this episodic mass loss oc¬ 
curring such a short time before core collapse suggests 
that it is connected to some instability in the latest 
stages of nuclear burning — most likely during Ne, O, 
and Si bur ning, but perhaps also C burning on longer 
timescales dOua taert fc Shiodd I2OI2I : IShiode fc Quataero 
l2014l : [Smith fc Arnettll2014l j. Some SNe have associted pre- 


(iFolev et al.ll2007l: Pastorello et al.l 

20071: ISmith et al.ll2010al: 

Mauerhan et al.ll2013al: ICorsi et al 

2 OI 4 I: iFraser et al.l 20131: 

Ofek et al. 2013a. 2014ah. Since stars mav alter their struc- 


ture shortly before core collapse, this makes it challenging 
to connect SNe to the properties of traditional populations 
of evolved stars. 

Eruptions driven by th e pulsational pair-instability 
(PPp lRakavv fc Shaviviri967l j have been suggested as a pos¬ 
sible phy sical cause of some o f this violent late-phase mass 
loss (e.g.- IWooslev et al]l2007ll . While the PPI could poten¬ 
tially be important in the most extreme and rare events 
like superlu minous SNe Iln that requ i re ejection of 10- 
20M(7) (e.g.. ISini th et al. ll2008al. l2010bl : lOfek et ^ l2014bl : 
iGal-YarJ l2012i r Zhang et al.l 2012|j . the requirements on 
the very high init ial mass of a PPI progenitor 100 M©; 
iHeger et al.ll200a l as well as a preference for low metallic- 
ity indicate that such events must be far too rare to ac¬ 
count f or the general populat i on of SNe Iln (| Smith et al.l 
2011^; ISmith fc Arnett! l2014l : ISmith. Mauerhan. fc Prietol 
2 OI 4 I). Multidimensional numerical simulations of the latest 


phases of nuclear burning in massive stars reveal instabili¬ 
ties that may be more widespread (i.e., in lower-mass stars 
of 20-30 M©), and might l ead to violent erupt i ve ma ss loss 
dArnett fc Meakinl 1201 il l. lOuataert &: Shiod3 (l2012l l sug¬ 
gested that these latest phases of Ne and O burning might 
induce wave-driven mass loss that could contribute (see also 
IShiode fc Quataerd 20lJ). A number of other possibilities 
may exist as well: Smith fc Arnett! (120141 ') pointed out that 
even if waves propagating out from the core are insufficient 
to cause hydrodynamic mass loss on their own, their dis¬ 
sipation in the envelope may cause a pulsation or inflation 
of the star, perhaps triggering a collision with a compan¬ 
ion star in a binary system that had previously been safely 
separated. Moreover, violent convection and mixing in the 
final phases may be able to trigger explosive nuclear burn¬ 
ing that could cause sudden energy deposition. In any case, 
the empirical fact that some SNe experience violent erup¬ 
tive mass loss in their latest burning phases indicates that 
the stellar structure may be significantly modified compared 
to the endpoints of one-di mensional stellar-evolution mod¬ 
els; ISmith fc Arnetj (120141 ') proposed that this could be an 
extremely important consideration in understanding core 


collapse (i.e., pre-explosion disruption of the core structure 
might make some stars easier to explode). 

SNe Iln contribute a fraction of 8-9% of all core-collapse 
SNe (ISmith et al.ll20lTal l in a volume-limited sample of large 
star-forming galaxies, and there are hints that they are as 
common or possibly more common at lower metallicity (see 
ISmith|[2014l l. It is interesting that ~ 10% of core-collapse 
SNe exhibit vivid warnings of their impending core collapse 
in the form of violent fits of dynamic mass loss, while the 
vast majority of normal SNe do not. An open question is 
whether explosions classified as SNe Iln are the only mas¬ 
sive stars that experience pre-SN instability, or if they are 
just the most extreme expression of a more generic instabil¬ 
ity operating in a larger fraction of massive stars. Whereas 
the PPI is indeed limited to a small fraction of the most 
massive stars, wave-driven mass loss, inflation of the star 
that triggers collisions with a companion star, and explosive 
burning may be more generic and may operate over a wide 
range of mass. Quantitative expectations of the energy de¬ 
position and mass loss from theory are still uncertain, but 
one can imagine that there is a substantial range in both. If 
there is a range in energy deposition from pre-SN instabil¬ 
ity, there may also be a continuum in CSM density, from the 
most violent SNe Iln down to normal winds. Anecdotally, it 
is already clear that SNe Iln exhibit a wide range of CSM 
mass with diverse radial density structure and geometry. 
They can eject as much as 10-20 M© in the cases of rare su- 
perlumin ous SNe Iln (see above), 0.1-1 M© in more common 
SNe Iln dKiewe et al.ll201 lOfek et ^l2013bl : iTaddia et al.l 
I2OI3II . and they even have CSM that could ar ise from the 
dense wi nds of extreme red superg i ants (RSCs; ISmith et al.l 
I2009al f9: iMauerhan fc Smithll2012 : IStritzinger et ah 2012| j. 


In this paper, we report evidence that the Type Iln phe¬ 
nomenon may extend to even lower mass-loss scales, in SNe 
that are recognised as Type Iln only if they are caught early 
enough after explosion, and then transition into more nor¬ 
mal SN types as they age. Thus, some other “normal” SNe 
might also experience pre-SN instability akin to SNe Iln, but 
less extreme in scale or more limited in duration. In the title 
of this paper and throughout the text, we refer to “normal” 
SNe as those which are not seen to be strongly interacting 
based on t heir visible-wavel ength spectra (SNe Ibc, Ilb, II- 
P, II-L; see lFiliDpen^ll997l for a review), where the visual- 
wavelength spectrum is dominated by the receding photo¬ 
sphere that is located within the fast freely expanding SN 
ejecta, and not in the shock/CSM interaction region. Thus, 
normal SNe in this context exclude Types Iln and Ibn. Signs 
of weaker CSM i nteraction have often been s e en in X-ray and 
radio emission (I Chevali er fc FranssorJ[l994l : iFransson et al.l 
ll996l:lMurase et al.ll2014fb 

We present visual-wavelength photometry and spec¬ 
troscopy of PTFlliqb taken from the time of first detec¬ 
tion up to ~ 1100 days afterward. PTFlliqb was discov¬ 
ered (iParrent et al.ll201ll ) at R = 16 .8 mag by the Palomar 
Transient Factory (PTF; iLaw et H] l2009l : lRau et al.l l2009h 
on 2011 July 23 (UTC dates are used throughout this pa¬ 
per) in an inner spiral arm of the nearby barred spiral galaxy 
NCC 151 (see Figure[T]), and was initially classified as having 
a Type Iln spectrum (as we discuss below, however, the spec¬ 
tral evolution became complicated). Based on the redshift 
of z = 0.0125 (3750 km s“^), we adopt a distance of 50.4 
Mpc (m — M = 33.51 mag), and a Milky Way line-of-sight 
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Figure 1. B-band image of NGC 151 taken on 2011 August 5 
with the Mount Laguna Observatory 40-inch reflecting telescope. 
PTFlliqb is indicated by the arrow. The local standards listed 
in Table [3] are marked. 


reddening of E{B — V)= 0.0284 mag (ISchlegel et al.l[l9^ . 
PTFlliqb resides in the bright inner region of NGC 151, in 
the ring surrounding the central bar, at a projected separa¬ 
tion of ~ 30” 7-8 kpc) from the host galaxy’s nucleus. 

Early reports indicated n ondetections of PTFlliqb in X - 
rays dOuimby et al.l 1201 il l and radio dHoresh et al.l I 2 OIIII . 
Later, PTFlliqb exhibited a marginal X-ray detection at 
~ 24 days after peak, and then another non-detection 28 
days after peak, which was suggested to indicate a relatively 
low ratio of the X-ray luminosity to visu al-wavelength lum i- 
nosity, as compared to other SNe Iln dOfek et al.ll2013al . 
Our new observations are presented in §2, and the light 
curve and spectral evolution are analyzed in §3. We com¬ 
pare PTFlliqb with SN 1998S in §4, and in §5 we present a 
simple CSM interaction model. In §6 we discuss the late-time 
Ha and in §7 we discuss the WR features seen in the early- 
time spectrum. Section 8 presents an overview and summary, 
and we discuss PTFlliqb in context with other SNe Iln and 
normal (noninteracting) Type II (II-P and II-L) SNe. 


Table 1. P48 photometry" 


MJD 

R (mag) 

cr (mag) 

55765.4 

16.862 

0.036 

55765.4 

16.834 

0.039 

55766.4 

16.135 

0.049 

55766.4 

16.073 

0.050 

55768.4 

15.608 

0.027 

55768.4 

15.603 

0.050 

55769.4 

15.462 

0.047 

55769.4 

15.458 

0.050 

55770.4 

15.370 

0.040 

55770.4 

15.323 

0.057 

55771.4 

15.291 

0.047 

55771.4 

15.265 

0.057 

55774.3 

15.218 

0.048 

55774.4 

15.194 

0.054 

55775.5 

15.147 

0.072 

55779.4 

15.194 

0.050 

55779.5 

15.178 

0.056 

55780.4 

15.169 

0.054 

55780.5 

15.180 

0.069 

55782.4 

15.204 

0.060 

55782.4 

15.193 

0.076 

55783.3 

15.247 

0.025 

55783.4 

15.258 

0.034 

55784.3 

15.274 

0.028 

55784.4 

15.274 

0.051 

55793.3 

15.482 

0.041 

55794.3 

15.486 

0.043 

55811.3 

15.789 

0.047 

55811.3 

15.776 

0.056 

55812.2 

15.779 

0.044 

55813.3 

15.822 

0.049 

55813.4 

15.807 

0.055 

55821.2 

15.909 

0.052 

55821.3 

15.899 

0.058 

55822.2 

15.921 

0.046 

55822.3 

15.923 

0.062 

55823.2 

15.921 

0.049 

55823.3 

15.908 

0.062 

55824.3 

15.939 

0.055 

55824.3 

15.897 

0.070 

55838.2 

16.054 

0.050 

55842.3 

16.082 

0.039 

55842.3 

16.085 

0.033 

55850.2 

16.137 

0.051 

55850.2 

16.134 

0.055 

55851.2 

16.154 

0.053 

55851.2 

16.160 

0.049 

55852.2 

16.144 

0.067 

55852.3 

16.159 

0.072 

55853.2 

16.169 

0.066 

55853.3 

16.186 

0.059 

55854.2 

16.131 

0.064 

55854.3 

16.137 

0.072 


“Mould R filter 


2 OBSERVATIONS 

2.1 Palomar 48-inch Discovery and Photometry 

We obtained R-band images of the PTFlliqb field includ¬ 
ing NGC 151 on 2011 July 22.37 with the Palomar 48-inch 
telescop e (P48) equipped w ith the refurbished CFHT12k 
camera (iRahmer et al.ll2008l L Subtraction of a stacked ref¬ 
erence image of the field with HOTPANT^ revealed a 
new transient source at coordinates a = 00*'34“'04’;84, 
S = —09°42”17'.'9 (J2000.0), with an astrometric uncer¬ 
tainty (relative to the USNO-Bl catalog: iMonet et al.ll200'^ l 
of ±150 mas in each coordinate. The transient was discov¬ 
ered 30 hours later by Oarical, an autonomous software 

^ http://www.astro.washington.edu/users/becker/hotpants.htm 


framework of the PTF collaboration (iBloom et ahllioiHl . It 
was classified correctly as a transient source (as opposed to 
a variable star), was further classified as a SN or nova, and 
was given the name PTFlliqb. No source was detected at 
this location with P48 in an image taken on 2011 Jul 17.5, 
to a 3cr limit of J? « 18.5 mag (see Figures [2] and [3]). 

Our photometric pi peline has been used in many pre¬ 
vious PTF papers (e.g . . [Firth et ah 1201,4 lOfek et ^1120121 . 
l2014al : |Pan et ^120141 : iLaher et al. 1201411 . and is based on 
image subtraction. We construct a deep reference image 
from data prior to the SN explosion, register this reference to 
each image containing the SN light, match the point spread 
functions (PSFs), perform the image subtraction, and then 
measure the SN flux using PSF photometry on the differ¬ 
ence images. The PSF is determined using isolated stars in 
the unsubtracted images, and the image subtraction uses a 
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nonp arametric pixelised kernel (similar to that in iBramicbl 
1200^ . The average PSF is then fit at the position of the SN 
event weighting each pixel according to Poisson statistics, 
yielding a SN flux and flux uncertainty. We flux calibrate 
our P48 lig ht curve to the Sloan Digital S ky Survey (SPSS ; 
lYork et al.fl200Cll l Data Release 10 (DR10: rAhn et Sll2014ll . 
The resulting magnitudes are given in Tableland are shown 
in Figures [5] and 13] 


2.2 Palomar 60-inch Photometry 


Upon discovery of PTFlliqb, the field was automatically 
inserted into the queue of the r obotic Palomar 60-inch tele¬ 
scope (P60; ICenko et al.l I 2 OO 6 II for multifilter (r' and i') 
follow-up observations. Images were processed using a cus¬ 
tom pipeline, and subtracted from reference images obtained 
in 2014 July using HOTPANTS. The resulting subtracted 
images were photometrically calibrated using nearby point 
sources from the SDSS DRIO dAhn et al.||2014[) ), so repo rted 
magnitudes are on the AB system ( Oke fc Gunn|[l983l l. A 
log of P60 photometry is provided in Table [2] The light 
curve, combining the P48 J?-band and the P60 r' and i' 
photometry, is shown in Figures [2] and [3| 

The general rate of fading in P60 ri data is very consis¬ 
tent with the P48 R-band data. However, there are some dif¬ 
ferences that can mostly be attributed to different amounts 
of Ho emission in the r (Sloan) and R (Mould) filters, and 
none in the i filter. We also include a number of i-band upper 
limits at late times around days 500-600. 


2.3 Mount Laguna BVRI Photometry 

All photometric data were collected using the Mount Laguna 
Observatory (MLO; Smith & Nelson 1969) 40-inch reflect¬ 
ing telescope, which is equipped with a 2048 x 2048 pixel 
CCD camera (manufactured by Fairchild Imaging Systems; 
CCD447) located at the //7.6 Cassegrain focus, providing 
a field of view of approximately 13.5' x 13.5' with 0.41" 
pixel” The “seeing,” estimated from the full width at half¬ 
maximum intensity (FWHM) of stars on the CCD frame, 
was generally ~ 2", and exposure times of 5 to 30 min¬ 
utes were typical for the observations, which were taken in 
the standard Johnson-Cousins (Johnson et al. 1966; Cousins 
1981) BVRI bandpasses. 

CCD frames were fiatfielded using either twilight-sky 
or dome fiats in the usual manner, and cosmic rays were re¬ 
moved using the L. A. Cosmic (van Dokkum 2001) algorithm. 
Considerable fringing remained in 7-band images that did 
not properly flatten, likely produced by the varying inten¬ 
sity of night-sky emission lines. This has minimal impact on 
the photometry since both the SN and the comparison stars 
were much brighter than the background. 

Figure [T] shows an MLO 73-band image of NGC 151 
taken on 2011 August 5, with 10 “local standards” identified 
in the field of PTFlliqb, which were used to measure the rel¬ 
ative SN brightness on nonphotometric nights. The absolute 
calibration of the field was accomplished on the photomet¬ 
ric night of 2012 November 19 by observing several fields 
of Landolt (1992) standards over a range of airmasses. The 
derived BVRI magnitudes of the stars are given in Tabled 
along with the estimated uncertainties. The transformation 


coefficients to the standard Johnson-Cousins systems were 
derived using the solutions from this night. We determined 
the instrumental magnitudes for the standards using aper¬ 
ture photometry with the IRAF0 DAOPHOT package (Stet¬ 
son 1987, 1991), which yielded colour terms for the MLO 
observations of the form 

B = 6-b0.093(H-U)-kC'B, 

V = V - 0.067 (B-V) + Cv, (1) 

R = r + 0.093 (V-R) + Cr, 

I = i + 0.005 iV-I) + Ci, 

where bvri are the instrumental and BVRI the standard 
Johnson-Cousins magnitudes. The terms Cb, Cv, Cr, and 
Cl are the differences between the zero-points of the instru¬ 
mental and standard magnitudes, determined for each ob¬ 
servation by measuring the offset between the instrumental 
and standard magnitudes and colours of local standards. 

We determined the instrumental magnitudes for 
PTFlliqb and the local standards using the point-spread 
function (PSF) fitting technique within DAOPHOT (Stet¬ 
son 1991, and references therein). We used only the inner 
core of PTFlliqb and the local standards to fit the PSF in 
order to reduce errors that can be introduced when there 
is a strong gradient in the background (e.g., Schmidt et al. 
1993). In practice, this core was generally set to be about 
the FWHM of a given image. While the fitting radius of 
the SN and comparison stars was varied from night to night 
to match the seeing, the sky background of the SN and lo¬ 
cal standards was always set to an annulus with a radius 
of 20-30 pixels (8'.'2-12'.'3) to maintain consistency through¬ 
out the observations. We next subtracted the mode of the 
sky backgrounc0 to derive the instrumental magnitudes for 
the SN and local standards. The transformation to the stan¬ 
dard Johnson-Cousins system was then accomplished using 
the colour terms listed in Equation 1 and the standard mag¬ 
nitudes of the local standards given in Tabled On the two 
nights where images in all four filters were not obtained, 
for the purpose of estimating the appropriate colour term 
the colour of PTFlliqb was estimated through interpola¬ 
tion or extrapolation from temporally nearby values; since 
the colour terms are small, this approximation should have 
little impact on the resulting photometry. 

The final reported photometry was accomplished by 
taking the simple mean of the values obtained using each 
of the available calibrator stars on a given night; not all 
stars were available on all nights, owing either to field-of- 
view limitations or saturation. The results of our photomet¬ 
ric observations are given in Table |T] and shown in Figure [31 
The reported uncertainties come from two sources. First, 
there is the photometric uncertainty reported by the error- 
analysis package in DAOPHOT from the statistics of the 
SN and background region. Second, there is uncertainty in 
the transformation to the standard system. We estimate the 


^ IRAF is distributed by the National Optical Astronomy Obser¬ 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

® See Da Costa (1992) for a discussion of the advantages of us¬ 
ing the mode of the background region rather than the mean or 
median. 
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Table 2. P60 r and i photometry of PTFlliqb® 


MJD 

Filt. 

mag 

cr 

MJD 

Filt. 

mag 

CT 


55766.449 

r 

16.134 

0.004 

55768.358 


15.948 

0.076 

55770.355 

r 

15.347 

0.004 





55774.344 

r 

15.230 

0.004 

55774.343 


15.445 

0.035 

55776.338 

r 

15.236 

0.004 

55776.337 


15.460 

0.069 

55776.340 

r 

15.217 

0.004 

55776.339 


15.464 

0.053 

55778.331 

r 

15.374 

0.004 

55778.330 


15.486 

0.107 

55778.334 

r 

15.255 

0.004 

55778.333 


15.441 

0.067 

55780.330 

r 

15.396 

0.003 





55782.321 

r 

15.357 

0.004 

55782.320 


15.460 

0.063 

55782.323 

r 

15.306 

0.004 

55782.322 


15.473 

0.097 

55792.482 

r 

15.417 

0.004 

55792.481 


15.720 

0.097 

55794.300 

r 

15.557 

0.004 





55794.300 

r 

15.704 

0.004 





55796.282 

r 

15.633 

0.004 

55796.281 


15.763 

0.049 

55798.300 

r 

15.839 

0.004 





55800.373 

r 

15.593 

0.004 

55800.372 


15.825 

0.046 

55801.400 

r 

15.779 

0.005 





55802.300 

r 

15.616 

0.004 





55803.264 

r 

15.821 

0.004 

55803.263 


15.948 

0.040 

55804.261 

r 

15.936 

0.004 

55804.260 


15.932 

0.070 

55805.259 

r 

15.815 

0.004 

55805.258 


15.986 

0.084 

55806.256 

r 

15.842 

0.005 






r 



55807.254 


15.935 

0.049 

55808.251 

r 

15.888 

0.054 

55808.250 


16.033 

0.066 

55809.400 

r 

15.807 

0.005 





55811.244 

r 

15.904 

0.005 

55811.243 


16.005 

0.058 

55812.239 

r 

15.965 

0.005 

55812.238 


16.090 

0.076 

55813.200 

r 

15.881 

0.006 





55814.233 

r 

15.725 

0.007 

55814.232 


16.151 

0.042 

55820.257 

r 

16.005 

0.006 

55820.256 


16.253 

0.059 

55821.216 

r 

15.887 

0.006 

55821.215 


16.248 

0.066 

55827.229 

r 

16.107 

0.005 

55827.228 


16.344 

0.061 

55828.223 

r 

16.060 

0.005 

55828.222 


16.379 

0.052 

55829.200 

r 

16.036 

0.006 





55830.200 

r 

16.265 

0.006 





55831.188 

r 

16.278 

0.158 

55831.187 


16.381 

0.050 

55832.211 

r 

16.165 

0.006 

55832.210 


16.341 

0.063 

55833.200 

r 

16.080 

0.006 

55833.181 


16.326 

0.048 

55834.179 

r 

16.1220 

0.006 

55834.178 


16.380 

0.086 

55835.177 

r 

16.076 

0.005 

55835.176 


16.358 

0.048 

55837.300 

r 

16.425 

0.005 





55838.168 

r 

16.551 

0.007 

55838.167 


16.616 

0.050 


r 



55847.144 


16.598 

0.058 

55848.141 

r 

16.310 

0.008 

55848.140 


16.524 

0.047 

55849.137 

r 

16.219 

0.008 

55849.136 


16.594 

0.062 

55850.139 

r 

16.210 

0.006 

55850.138 


16.639 

0.058 

55851.133 

r 

16.135 

0.006 

55851.132 


16.597 

0.085 

55852.129 

r 

16.176 

0.006 

55852.128 


16.630 

0.064 

55852.200 

r 

16.208 

0.005 





55853.130 

r 

16.323 

0.006 

55853.129 


16.640 

0.067 

55854.129 

r 

16.330 

0.006 

55854.128 


16.687 

0.086 

55856.100 

r 

16.324 

0.006 





55857.100 

r 

16.269 

0.006 





55858.100 

r 

16.154 

0.006 






r 



55859.123 


16.813 

0.053 

55862.113 

r 

16.555 

0.007 

55862.112 


16.852 

0.062 

55863.111 

r 

16.481 

0.007 

55863.110 


16.792 

0.041 

55864.100 

r 

16.401 

0.007 





55866.110 

r 

16.422 

0.007 

55866.109 


16.868 

0.069 

55867.100 

r 

16.467 

0.009 





55868.105 

r 

16.693 

0.008 

55868.104 


16.952 

0.060 

55869.100 

r 

16.509 

0.012 





55874.297 

r 

16.533 

0.014 

55874.296 


17.083 

0.065 

55884.100 

r 

16.892 

0.012 





55889.112 

r 

17.194 

0.010 

55889.111 


18.203 

0.069 

55891.296 

r 

17.411 

0.016 





55893.157 

r 

17.576 

0.016 

55893.156 


18.604 

0.075 


r 



55895.094 


18.620 

0.074 

55902.247 

r 

17.524 

0.024 

55902.245 


18.907 

0.081 

55903.112 

r 

17.927 

0.017 

55903.110 


18.836 

0.058 


r 



55904.109 


18.903 

0.078 


r 



55905.105 


(19.76) 


55906.276 

r 

17.595 

0.020 

55906.274 


19.040 

0.095 

55916.133 

r 

17.537 

0.016 

55916.131 


19.120 

0.099 


r 



56290.090 


(21.69) 


56296.100 

r 

20.647 

0.332 





56297.140 

r 

20.617 

0.131 

56297.139 


(21.86) 


56298.115 

r 

20.723 

0.109 

56298.113 


(22.33) 


56301.084 

r 

20.673 

0.111 

56301.081 


(21.91) 


56302.107 

r 

21.423 

1.151 

56302.106 


(19.96) 


56324.100 

r 

20.660 

0.211 

56324.128 


(22.19) 
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r 



56325.129 


(20.21) 


r 

21.129 

0.165 

56328.113 


(21.92) 


56330.112 

r 

21.142 

0.173 

56330.110 


(22.03) 



Magnitudes in narenthesis indicate nnner limits. 
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Table 3. MLO Magnitudes of Local Standards 


Star 

B (o-b) 

V (a-y) 

R (o'h) 

I (<^i) 

1 

13.822(0.040) 

13.258(0.022) 

12.873(0.033) 

12.563(0.033) 

2 

17.192(0.052) 

16.319(0.027) 

15.709(0.039) 

15.192(0.040) 

3 

15.522(0.040) 

14.898(0.022) 

14.408(0.033) 

14.019(0.033) 

4 

16.091(0.040) 

14.796(0.022) 

13.859(0.033) 

13.129(0.033) 

5 

16.856(0.045) 

15.714(0.023) 

14.949(0.035) 

14.250(0.034) 

6 

15.253(0.040) 

14.511(0.022) 

14.025(0.033) 

13.591(0.033) 

7 

14.883(0.040) 

14.214(0.022) 

13.823(0.033) 

13.446(0.033) 

8 

17.893(0.065) 

16.660(0.029) 

15.858(0.042) 

15.055(0.042) 

9 

14.702(0.040) 

14.097(0.022) 

13.734(0.033) 

13.339(0.033) 

10 

11.928(0.040) 

11.416(0.022) 

11.093(0.033) 

10.778(0.033) 


Table 4. MLO Photometric Observations of SN PTFlliqb 


UT Date“ 

Dayi> 

MJD 

B (ctb) 

V (try) 

R (o'h) 

t (o-j) 

2011-08-03 

11.06 

55776.46 



15.049(0.019) 

15.031(0.014) 

2011-08-06 

14.07 

55779.47 

15.334(0.064) 

15.260(0.025) 

15.051(0.082) 

14.949(0.034) 

2011-09-01 

39.99 

55805.39 



15.619(0.064) 


2011-09-17 

56.02 

55821.42 

17.036(0.045) 

16.447(0.097) 

15.902(0.036) 

15.623(0.039) 

2011-09-23 

61.95 

55827.35 

17.159(0.030) 

16.546(0.052) 

15.998(0.040) 

15.748(0.025) 

2011-09-28 

67.00 

55832.41 

17.256(0.069) 

16.518(0.042) 

15.976(0.035) 

15.727(0.033) 

2011-10-19 

87.91 

55853.31 

17.582(0.037) 

16.877(0.063) 

16.160(0.055) 

15.985(0.033) 

2011-10-27 

95.88 

55861.28 

17.738(0.075) 

16.933(0.063) 

16.258(0.037) 

16.098(0.041) 

2011-11-02 

101.86 

55867.26 

17.790(0.058) 

17.073(0.041) 

16.369(0.046) 

16.190(0.027) 

2011-11-15 

114.89 

55880.29 

18.130(0.041) 

17.525(0.068) 

16.717(0.024) 

16.546(0.027) 

2011-11-18 

117.88 

55883.28 

18.191(0.043) 

17.644(0.036) 

16.832(0.032) 

16.710(0.027) 

2011-11-30 

129.83 

55895.23 

18.599(0.034) 

18.485(0.023) 

17.533(0.026) 

17.719(0.031) 


Note: All photometric observations were made with the Mount Laguna Observatory 40-inch telescope. 

“yyyy-mm-dd. 

*’Days since first detection, 2011-07-22.9 UT (MJD 2,455,765.4). 


transformation error by taking the standard deviation of 
the spread of the standard magnitudes obtained using each 
of the local standard stars. The photometric and transfor¬ 
mation errors were then added in quadrature to obtain the 
uncertainty reported in Table |4] in nearly every case, the 
total error was dominated by the uncertainty in the trans¬ 
formation. 

2.4 Swift/UVOT Ultraviolet Photometry 

The Ultra-Violet Optical Telescope (UVOT; iRoming et al.l 
|2005^ onboard Swift observed PTFlliqb in the UVWl, 
UVM2, and UVW2 filters beginning on 2011 July 24. We 
photometered the UVOT data using standard procedures 
within the HEASoftQ environment. We used a 3” diameter 
aperture to extract flux from the transient, and subtracted 
coincidence-loss-corrected count rates from un derlying host 
emiss ion obtained from images in 2014 (e.g., iBrown et al.l 
I2OO9I ). Phot ometric cali b ration was performed using the 
recipes from IPoole et al] (l2008l ). The resulting magnitudes, 
all corrected to the AB system, are displayed in Table [S] 

2.5 Spectroscopy 

PTFlliqb was brightening quickly at the time of discovery 
and was presumably caught very early after explosion (see 
below), so we quickly initiated spectroscopic follow-up 
observations starting on day 2 after first detection (1 day 

^ http://heasarc.nasa.gov/lheasoft/ 


after discovery). We obtained several epochs of optical 
spectroscopy of PTFlliqb, which were densely sampled 
at early times (see Fig. [21 where epochs of spectroscopic 
observations are plotted above the light curve). We ob¬ 
tained spectra of PTFlliqb using a number of different 
facilities, including the Bluechannel (BC) spectrograph on 
the 6.5-m Multiple Mirror Telescope ( MMT), the Multi- 
Ob je ct Double Spectrograph (MODS; iBvard &: O’BrienI 
2OOOI) on the LBT. the Double Beam Spectrograph (DBSP; 
Oke fc GunnI 1 19821 ) on the Palomar 200-inch telescope 
(P200), the Intermediate dispersion Spectrograph and 
Imaging System (ISI^I) on the 4.2-m William Herschel 
Telescope (W HT), CMOS o n Gemini South, the Kast 
spectrograph jMiller fc Stonelll993l ) on the Lick 3-m Shane 
reflector, the Lo w-Resolution Imaging Spectrometer (LRIS; 
lOke et ^Il995l) mounted on the 10-m Keck I telescope, and 
the Deep Irnaging Multi-Object Spectrograph (DEIMOS; 
[Faber et al.l l2003li on Keck 11. Details of the spectral 
observations are summarised in Table [5] and the spectra 
will be released on the Weizma nn Interactive Supernova 
Data Repository (WISeREP; lYaron fc Gal-Yam I2OI2I : 
http://www.weizmann.ac.il/astrophysics/wiserep/). 

The slit w as _ ^ways oriented at the parallactic angle 

llFilippenkdll983) . and the long-slit spectra were reduced 
using standard procedures. Final spectra are shown in 
Figure[4l although a few epochs of spectra are not displayed 
in this figure because they were taken close in time to 
another spectrum that appears very similar. The earlier 

® http://www.ing.iac.es/PR/wht_info/whtisis.html 
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Figure 2. Top: Apparent magnitudes for the early part of the light curve for PTFlliqb, with P48 i2-band in solid red dots, and with P60 
r' and i' bands shown with unfilled circles in red-orange and black, respectively (see Tables^ and [2jl- Swift/UVOT ultraviolet photometry 
in the UWl, UM2, and UW2 filters is shown in magenta, blue, and purple, respectively (see Table igi. UVOT magnitudes are offset by 
-\-2 for clarity. BVRl photometry from Mount Laguna Observatory (MLO) is shown in blue, green, red, and black squares (Table [4|), 
with all four MLO bands offset by -h2 mag for clarity. These are all observed apparent magnitudes (i.e., not corrected for Milky Way 
reddening). Bottom: Colour evolution of PTFlliqb in the UV (UVW1-UV W2; purple) and visib le wavelengths (B — V; black squares). 
For comparison, B — V colours for the SN II-P SN 1999em ft hick grev line: [Leonard et al.|[20Q3i . SN II-L 2003hf (orange circles; from 
iFaran et al.||2of4) . and the SN Iln SN 1998S (green diamonds: iFassia et al.ll2000li are shown. 


Table 5. Swift/UVOT photometry of PTFlliqb 


MJD 

UM2 (mag) 

a- (mag) 

UWl (mag) 

<7 (mag) 

UW2 (mag) 

cr (mag) 

55766.1 

15.14 

0.05 

15.28 

0.05 

15.15 

0.06 

55768.0 

14.87 

0.05 

14.91 

0.05 

14.93 

0.06 

55770.9 

15.08 

0.05 



15.33 

0.06 

55773.7 

15.57 

0.05 

15.39 

0.05 

15.87 

0.06 

55775.6 

16.03 

0.05 

15.70 

0.05 

16.37 

0.06 

55777.2 

16.39 

0.05 

16.01 

0.05 

16.77 

0.06 

55781.9 

17.19 

0.07 

16.70 

0.05 

17.55 

0.08 

55787.1 

18.19 

0.07 

17.58 

0.06 

18.51 

0.09 

55797.1 

19.07 

0.11 

18.30 

0.07 

19.44 

0.13 

55803.4 

19.28 

0.12 

18.59 

0.09 

19.87 

0.18 

55807.1 

19.54 

0.13 

18.80 

0.10 

19.87 

0.16 
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0 50 100 150 200 300 400 500 600 700 800 

Days after first detection 

Figure 3. Apparent and absolute magnitude light curve for PTFlliqb, with P48 R in solid red dots, and with P60 r' and shown 
with unfilled circles in red-orange and black, respectively (see Tables^ and [3. The arrows at late times are P60 upper limits, whereas 
the solid rectangles at late times (after day ~ 150) represent pseudo-photometry estimated from spectra as described in the text. All 
photometry for PTFlliqb in the three filters has been corrected for Milky Way reddening of 0.070, 0.073, and 0.054 mag, for R, r', and 
z', respectively. Epochs when we obtained visual-wavelength spectra are noted with co loured hash-marks a t the top (see Table 3 and 
Figure[^. For refer ence, we also show J?-band light curv es of the normal SN II-P 1999em l|Leonard et al.l2Q02li . the unusual SN Iln 2010jp 
llSmith et al.ll2012t) . and SN 1998S llFassia et al.|[200oh . scaled as they would appear if they were at the same distance as PTFlliqb. 
The dashed line is representative of the ®®Co decay rate, matched to SN 1999em. Note that the time axis is not uniform; we include a 
change in scale at 200 days to better facilitate the display of the early-time plateau phase on the same plot as the more sparsely sampled 
late-time decay. 


epochs in Figure|3]have a blackbody plotted in grey for each 
observed spectrum; these correspond to the temperatures 
listed as Tbb in Table |6l although these are intended only 
as a rough relative comparison because they depend on 
the adopted reddening. Nevertheless, inferred temperatures 
around 7000 K in the plateau phase suggest that these 
quoted temperatures are not wildly in error. Details of the 
Ha line profile are shown in Figure [5] 

Since our P48 and POO photometry covers early times 
well, but includes poor sampling after day 150 (and only up¬ 
per limits in the i' band on days 500-600), we supplement 
our imaging photometry by using our flux-calibrated spec¬ 
tra. In spectra taken after day 150, we measured integrated 
fluxes within specific wavelength ranges to simulate the flux 
that would be observed in the R (5700-7300 A), r' (5500- 
6700 A), and i' (6900-8100 A) filters, and then converted 
these fluxes to magnitudes and added them to Figure [3] as 
solid rectangles with the same symbol colours as the imag¬ 
ing photometry. The uncertainty here is difficult to quantify. 


since the main source of error is the exact positioning of the 
target and associated standard star within the slit, but we 
adopt representative errors of ±0.4 mag in Figure |3] (at ear¬ 
lier times when our photometry overlaps with spectra, the 
spectrophotometry generally agree to 0.4 magor better). As 
discussed in more detail below, it is important to recognise 
that at the redshift of PTFlliqb, the main difference be¬ 
tween the R and r' filters is that R includes the full Ha 
emission-line profile, whereas the red end of the r' filter cuts 
off the red wing of PTFlliqb’s Ha line. It is also important 
to note that the pseudo-photometry from spectra subtracts 
the background by sampling host-galaxy light on either side 
of the SN, whereas the imaging photometry utilises back¬ 
ground subtraction with a SN-free reference image. Thus, if 
there is an underlying star cluster or other coincident source, 
it will be included in the spectral photometry but will be ab¬ 
sent in the imaging photometry. This may be important at 
late times when PTFlliqb is faint, so the true brightness 
of PTFlliqb may be below the level indicated by the spec- 
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Table 6. Spectroscopic Observations of PTFlliqb 


Date 

Tel./Instr 

Day 

AA(A) 

Tbb(k) 

2011Jul24 

GS/GMOS 

2 

3500-7400 

25,000 

2011Jul24 

P200/DBSP 

2 

3520-9825 

25,000 

2011Jul27 

WHT/ISIS 

5 

3500-9500 

16,000 

2011Jul30 

WHT/ISIS 

8 

3500-9500 


2011 AugOl 

Keck2/DEIMOS 

10 

5150-7800 


2011 Aug 04 

Lick/Kast 

13 

3436-9920 

10,500 

2011 Aug 05 

Lick/Kast 

14 

3436-9920 

10,500 

2011 Aug 06 

P200/DBSP 

15 

3520-9825 

10,000 

2011 Aug 07 

WHT/ISIS 

16 

3500-9500 

10,000 

2011 Aug 28 

P200/DBSP 

37 

3520-9825 

6,500 

2011 Sep 21 

P200/DBSP 

61 

3520-9825 

6,500 

2011 Sep 26 

Keckl/LRIS 

66 

3200-7350 


2011 Sep 29 

Keck2/DEIMOS 

69 

4500-9300 

6,500 

2011 Oct 25 

Lick/Kast 

95 

3436-9920 

5,600 

2011 Nov 26 

Keckl/LRIS 

127 

3340-10300 

7,500 

2011 Nov 26 

Lick/Kast 

127 

3436-9920 

7,500 

2011 Dec 02 

Keckl/LRIS 

134 

3300-7400 

7,500 

2011 Dec 18 

Lick/Kast 

150 

3436-9920 

7,500 

2011 Dec27 

Lick/Kast 

159 

3436-9920 

7,500 

2012 Jan 02 

MMT/BC 

165 

5550-7500 


2012 Jan 19 

MMT/BC 

182 

3820-9000 

7,500 

2012 Jan 26 

Bok/SPOL 

189 

3900-7600 

7,500 

2012 Jan 30 

MMT/BC 

193 

5550-7500 


2012 Jun 16 

Keckl/LRIS 

331 

3330-10100 


2012 Jun 28 

MMT/BC 

343 

5550-7500 


2012 Dec 18 

Keckl/LRIS 

516 

3400-10200 


2013 Jan 14 

MMT/BC 

543 

5550-7500 


2013 Jul 11 

Keck2/DEIMOS 

721 

4400-9600 


2013 Sep 04 

MMT/BC 

775 

nondetect. 


2013 Sep 29 

LBT/MODS 

801 

5500-10000 


2014 Jul 29 

Keck/LRIS 

1104 

3400-10200 



trophotometry. Nevertheless, these measurements provide a 
useful guide for the late-time evolution. 


3 RESULTS 
3.1 Light Curve 

The multiband light curves and colour evolution for the early 
phases of PTFlliqb are shown in Figure [21 while Figure [3] 
shows a subset of the data on an absolute-magnitude scale 
and extended to later times, as compared to a few previously 
observed SNe. The light curve displays a rapid rise in only 
~ 10 days to a peak luminosity of about —18.4 mag. This 
is followed by a decline that flattens out in a plateau or 
shoulder at about —17 mag, dropping oft the plateau after 
100-120 days. 

To the extent that flux proportional to is a valid 
description of the early rise, the dotted curve in Figure O 
suggests that explosion occurred at roughly —3.3 days (i.e., 
3.3 days before day 0, defined here as the first detection on 
2011 July 22.37). This is consistent with the upper limit 
of 18.5 mag (unfiltered) for our nondetection on 2011 July 
17.47 (Parrent et al. 2011), indicated by the red arrow in 
Figure O Since this upper limit corresponds to an absolute 
magnitude of around —15 mag, it does not place strong con¬ 
straints on the luminosity or duration of a possible pre-SN 
eruption. The relatively quick rise to peak in only a few 
days seems to suggest an extended progenitor star similar 
to a RSG, which is a topic we discuss in more detail later. 

The visible-wavelength colour evolution of PTFlliqb 
(see Figure mi starts out very blue and becomes steadily 
redder for the first 60-70 days, and then seems to level off 
a.t B — V ~ 0.75 mag before becoming more blue again at 



Figure 4. Low/moderate resolution spectra of PTFlliqb (see 
Table |6]l , with early times at the top and later times at the bot¬ 
tom. The four main phases discussed in the text are marked with 
changing colours: early times with CSM interaction (black), the 
plateau (blue), the early nebular phase (black), and the late nebu¬ 
lar phase dominated again by CSM interaction (magenta). Strong 
telluric features are marked with “®”. 

late times. PTFlliqb never gets as red as a normal SN II-P, 
like SN 1999em, probably owing to the continued influence 
of CSM interaction. The B — V colour evolution of PTFlliqb 
is quite similar to that of SN 1998S, although at early times 
it is also similar to a Type II-L event like SN 2003hf (the 
photometry here is from the recent study of SNe II-L by 
iFaran et al.ll201^ . We demonstrate later that the spectra 
around days 30-60 most closely resemble SNe II-L at simi¬ 
lar times. The UV colours are very blue at early times, but 
quickly redden as the UV luminosity plummets and the pho¬ 
tosphere cools rapidly in the hrst 30-40 days after explosion. 

There is little change in colour or luminosity during 
the latter part of the plateau phase from days ~ 50 to 
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120. After the drop from this plateau, PTFlliqb exhibits 
a nebular phase that (at least in the i' band) roughly fol¬ 
lows the luminosity decline expected for radioactive decay 
from ®®Co (see Figure [Sjl. While PTFlliqb is qualitatively 
similar to the light-curve behaviour of SNe Iln-P in this 
respect (iMauerhan et al]l2013hl : ISmithll2013bh . the drop af¬ 
ter the plateau was not nearly as severe, and the late-time 
radioactive-decay luminosity was not as low. Also, as we 
discuss below, during the decline from peak and the plateau 
phase, the spectrum does not show a strong Type Iln signa¬ 
ture (resembling SNe 1994W and 2011ht), but more closely 
resembles spectra of SNe ITP and II-L. Thus, its spec¬ 
tral evoluti on was unlike the rath er homogeneous class of 
SNe Iln-P (IMauerhan et al]|2013lJl . which exhibit a persis¬ 
tent SN Iln spectrum throughout their evolution. For these 
reasons we do not classify PTFlliqb as a Type Iln-P event, 
although we acknowledge that SN classification criteria may 
be inadequate if an object morphs through different types 
as it evolves. 

Based on an analysis of the spectra (below), it is clear 
that the i' band traces the true continuum luminosity much 
better than the R or r' bands. After the drop of the plateau 
in the light curve at 100-120 days, the red filters are dom¬ 
inated by strong Ha emission from CSM interaction, to 
a different degree in each filter because of their different 
wavelength responses (see above). The i' band does not 
include the bright Ha line. From the time of the plateau 
drop until about day 500, the i' band is consistent with ra¬ 
dioactive decay from ®®Co, and in fact has about the same 
luminosity and decline rate as the prototypical Type H- 
P SN 1999em. This would imply a synthesised ®®Ni mass 
of roughly 0.03 M|T) (for SN 1999e m, this was esti mated as 
0.02 M© bv lElmhamdi et al.ll2003l and 0.036 Mq bv lUtrobinI 
[20 o 3), and hence, a normal core-collapse SN H-P. 

As discussed later, we attribute the excess luminosity 
(compared to SN 1999em) in the first half of the plateau to 
CSM interaction, as well as the strong late-time Ha emis¬ 
sion. At the very latest phases (after day 500) the true un¬ 
derlying continuum level is very uncertain, since photometry 
from template-subtracted images only provides upper lim¬ 
its, whereas the pseudo-photometry estimated from spectra 
does not correct for a possible underlying star cluster or 
other faint coincident source, and may therefore be an over¬ 
estimate of PTFlliqb’s luminosity by an amount that is 
difficult to quantify. 


3.2 Spectral Evolution 

Here we provide a brief overview of the main stages in the 
spectral evolution of PTFlliqb. The spectra are discussed 
in more detail below in the comparison with SN 1998S, as 
well as in sections discussing the late-time Ha and early-time 
WR features specifically. We highlight three main stages: (1) 
at very early times around peak luminosity, (2) days 20-100, 
corresponding to the latter part of the plateau in the light 
curve, and (3) during the nebular phase after day ~ 120. 

3.2.1 Early Spectral Evolution 

For the first < 20 days after explosion, the spectrum of 
PTFlliqb was characterised by a smooth blue continuum 


matched by a blackbody that cools quickly from about 
>25,000 K on day 2 to about 10,000 K on day 16. The spec¬ 
trum on day 2 shows a narrow (unresolved at FWHM « 500 
km s“^ at this epoch) Ha line core with broader and sym¬ 
metric Lorentzian wings extending to about ±2000 km s“^, 
and also shows a comparably strong “Wolf-Rayet” emission 
feature (a combination of He ll A4686 and C lll/N ill). By 
day 5 these features have mostly faded; thereafter, only nar¬ 
row Balmer lines remain atop the smooth blue continuum. 
Similar WR features have been seen in a few other SNe at 
early times; this is discussed more in §7. 

3.2.2 Plateau Spectral Evolution 

Starting on day 16, we begin to see an underlying broad com¬ 
ponent of Ha that continues to grow in strength as the nar¬ 
row component fades (Figures |3] and [S^). This broad compo¬ 
nent arises in the underlying fast SN ejecta. There is some 
persistent very faint narrow Ha emission during the plateau 
phase, but it is difficult to ascertain how much of it may 
be contributed by underlying H ll regions, since its strength 
varies from one spectral epoch to the next (which could be 
caused by true intensity changes, or by various amounts of 
surrounding H ll region emission being included in the slit). 
In any case, the narrow component of Ha is faint (a few per¬ 
cent of the total Ha flux). On the other hand, the radius of 
the CSM interaction front in our model (see below) and the 
SN ejecta photospheric radius are comparable at this epoch, 
so perhaps the CSM interaction intensity remains steady, 
but much of the shock interaction region is overtaken by the 
opaque SN ejecta photosphere (see below). 

After about a month, the featureless hot blue contin¬ 
uum from the earlier phase transforms into a spectrum dom¬ 
inated by broad absorption and emission proHles, as in nor¬ 
mal ejecta-dominated SNe like SNe H-L and II-P. Besides 
the broad Ha line with a width of ~ 10,000 km s”*^ (FWZI), 
we see nu merous absorption f eatures in the blue typical of 
SNe H-P (jPessart et al.l[2013l i. as well as the broad Ca ll 
near-infrared (IR) triplet in emission. Some of the absorp¬ 
tion (especially Ha and Ca ii) appears muted compared to 
canonical SN II-P atmospheres. Indeed, PTFlliqb’s rela¬ 
tively weak P-Cygni absorption in Ha and Ca ll more closely 
resembles that of SNe II-L than SNe H-P, while at 4000- 
6000 A, SNe H-P and II-L look qua litatively very similar 
(see Figure [6t. ICutierrez et ah] ll2014h discuss the diversity 
in the ratio of Ha absorption to emission (a/e) in SNe H- 
P and ILL, and PTFlliqb resembles objects in that study 
with lower a/e vaues. So overall, even though the light curve 
shows a clear underlying plateau, the spectrum during the 
plateau phase looks more like that of SNe ILL; we suspect 
that this may be due to reionisation of the outer SN ejecta 
by CSM interaction. As we discuss below, the SN ILL-like 
spectral appearance is shared by SN 1998S at similar times 
in its evolution. 

3.2.3 Late-Time Spectral Evolution 

After the drop in luminosity at the end of the plateau (af¬ 
ter days 100-120), the Ha emission from PTFlliqb starts 
to become quite interesting, exhibiting an asymmetric and 
multipeaked profile (Figure[5]D). During this stage, Ha devel¬ 
ops a strong and somewhat broad blueshifted emission peak 
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Figure 5. Details of the Ha line profile of PTFlfiqb. (a) The left panel shows relatively early-time data, from the time of discovery 
until the beginning of the nebular phase. These spectra are normalised to the continuum level, and the time sequence of line strength 
(including broad components) leads to spectra that generally go up with time on this plot, (b) The right panel shows late-time spectra 
during the nebular phase from days 165 onward. Short vertical lines (dotted) mark velocities measured for relatively well-defined red and 
blue peaks (see text and Figure [8)l. 


centered at —3000 to —3500 km s“^. The red side of the line 
shows no similar peak, but instead begins to drop off sharply 
at around +3000 km s“^. This asymmetric blueshifted Ha 
profile changes little over the time period from days 120 to 
200 (Figure [5^). 

Blaming this strong blueshifted asymmetry of Ha on 
dust or opaque SN ejecta (blocking the receding side) would 
seem appealing, were it not for the fact that the asymmetric 
Ha profile then transitions from a blueshifted asymmetry to 
a much stronger redshifted peak between days 331 and 516 
(FigurelSlr). In the day 516 spectrum and afterward, the red 
peak of Ha remains stronger than the blue peak (even if dust 
that was previously blocking the red side were destroyed, it 
would not make the red side brighter; true asymmetry in 
the CSM density is needed). We note that a delay in light¬ 
crossing time from the back to the front of the SN cannot 
explain this shift from blue to red peaks, since the light¬ 
crossing time is less than 10 days at radii of ~ 500 AU. 

Interestingly, the strength of the blueshifted peak coin¬ 
cides temporally with a lingering blue continuum and broad 
SN ejecta lines like the Ca ll IR triplet. These are strong 
in the early phases of the nebular decay, as in typical SNe 
H-P. From day 516 onward, however, the blueshifted peak 
of Ha, the blue continuum, and the broad ejecta lines all 
fade in tandem. When the red peak is strongest (after day 
500), the overall spectrum is dominated by Ha with very lit¬ 
tle continuum contributing to the observed spectrum. This 
probably signifies that the spectrum at days 120-200 has 


a substantial contribution from the normal radioactivity- 
powered nebular emission from the inner ejecta, and this is 
when the i-band light curve and the Ha line luminosity (see 
below) both follow the ®®Co decay rate. At the latest times 
(after 300-500 days), the spectrum is dominated by CSM 
interaction because the radioactivity has faded much faster 
than the ongoing CSM interaction luminosity. 


A narrow emission component of Ha is seen at several 
epochs, most likely tracing the preshock CSM. Our highest- 
resolution spectra are those epochs obtained at late times 
with the MMT using the 1200 line mm“^ grating, which 
are most suitable to constrain the preshock wind speed. 
Examining epochs that appear to have less contamination 
from surrounding H ll regions (days 37, 61, 69, 95, 127, 150, 
165; probably owing to better seeing, evidenced by weaker 
[N ll] emission adjacent to Ha), we find that the narrow 
component of Ha has Gaussian FWHM values of ~ 80 km 
s“^, comparable to the spectral resolution of the data. The 
preshock wind is therefore much slower than that of a typical 
WR star (~ 10^ km s”*^ or more), or that of a blue super¬ 
giant or luminous blue variable (LBV; 100-500 km s“^), and 
clo ser to values seen in cool red or yellow supergiants (YSGs; 
see lSmi tbll20l4) . We quote 100 km s ^ as a fiducial value for 
the CSM expansion speed below, but the actual wind speed 
and mass-loss rates may be smaller by a factor of 2-3. 
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Figure 6. Comparison of visible-wavelength spectra of PTFlliqb 
on days 37 and 61 (blue) with those of SNe II-L (SN 2003hf; 
orange) and II-P (SN 1999em; red). For SN 1999em, we show the 
observed day 45 spectrum, as well as that same spectrum added to 
a blackbody to dilute the strength of absorption features. These 
spectra of SNe 2003hf and 1999em wer e obtained from th e UC 
Berkeley SN databas e, and appeared in iFaran et al. I ll2014l and 
iLeonard et alJ <l2002h . respectively. 


4 A COMPARISON WITH SN 1998S 


Readers familiar with SNe Iln may recognise that our de¬ 
scription of the spectral evolution of PTFlliqb in §3.2 
sounds remarkably similar to that of the classic well- 
observed Type Iln object SN 1998S. Figure [7] compares the 
spectral evolution of the two SNe, with visual-wavelength 
spectra interleaved in chronolog ical order after explo sion. 
The spectra of SN 1998S from [Leonard etay_([ 200(ji_axe 
taken from the UC Berkeley SN database ( Silverman et al.l 
|2012|~). and the othe rs are downloaded from the WISeREP 
l Yaron fc Gal-Ya^ 201 Jl d atab ase, originati n g fro m pa¬ 
pers by iFassia et al.l (2003) and IPozzo et akl (12004) . The 
days listed in Figure [7| are relative to discovery as noted 
by those authors; the explosion date of SN 1998S is not 
known precisely, but may be a few days before discovery 
([Leonard et al.li2000h . 

Basically, the spectra of the two SNe in Figure (3 are 
almost identical in their evolution, sufficient to claim that 
PTFlliqb is a near twin of SN 1998S. Both SNe proceed 
through all the same spectral changes outlined above at 
roughly similar times. We note two exceptions to this twin- 
hood. 

First, spectral signatures of CSM interaction are gen¬ 
erally stronger in SN 1998S at early times, consistent with 
its significantly higher peak luminosity (Figure O- In its 
day 2 spectrum, PTFlliqb shows an almost identical broad¬ 
winged Lorentzian Ha profile and Wolf-Rayet (WR) features 
to those that wer e first described in de tail for SN 1998S 
by [Leonard et al.l l[200Cl[f . [Chugaj ([200 ih noted how these 
can arise from CSM interaction in the opaque inner wind. 
An alternative explanation for the early-time spectrum may 



Figure 7. The spectral evolution of PTFlliqb compared to that 
of SN 1998S. All spectra of PTFlliqb are black, and all spectra of 
SN 1998S are red. Spectra of PTFlliqb for various epochs are the 
same as in Figure]^ Spectra of SN 1998S are fronjLeonardet_^ 
ll2000h (days 4, 25, 137, 31 2, 375, and 494'). [fissia et al.l ll200ll'l 
(days 12, 72, and 97), and IPozzo et all 112004^ (day 1093b The 
days for SN 1998S are as listed by those authors; this could cause 
offsets of a few days compared to PTFlliqb, which may be rele¬ 
vant at early times when comparing the two SNe, but which will 
not be significant at later times. 


be ionisation of a dense inner wind by a UV flash from 
shock breakout, as wa s hypothesised to occur in SN 2013cu 
(|Gal-Yam et al.l 120141) . It is difficult to distinguish between 
these two hypotheses with the available data, but in either 
case, PTFlliqb requires an opaque inner wind at a radius of 
~ 10 AU. An argument in favour of CSM interaction being 
important is that the wind density parameter required for 
the inner wind is very similar to that in the outer wind re¬ 
quired by the late-time CSM interaction when the emission 
is optically thin (see §7). In any case, the key difference as 
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compared to SN 1998S is that in PTFlliqb these features 
are weaker and they fade more qnickly. In the day 2 spec- 
trnm of PTFlliqb, we measure EW(Hq) = 13.5 ± 0.8 A 
and EW(WR) = 29.6 ± 1.0 A (emission-line EWs are posi¬ 
tive). The WR bnmp and strong Lorentzian wings of Ha are 
only seen in the day 2 spectrum and they fade completely by 
day 5 in PTFlliqb, leaving only a narrow Ha component. In 
SN 1998S, by contrast, these features are still quite strong i n 
the days 3, 4, and 5 spectra shown by [Leonard et all 1I2OOOII . 
In the day 4 spectrum of SN 1998S in Fignre[7l we measure 
EW(Ha) = 55T1.0 A and EW(WR) = 39T1.0 A. Note that 
because of the higher luminosity of SN 1998S, the line lumi¬ 
nosities are abont an order of magnitude larger in SN 1998S. 
It is also interesting to note that the EW(WR)/EW(Ha) ra¬ 
tio is higher in PTFlliqb, perhaps reflecting the fact that 
the spectrum was obtained sooner after explosion than for 
SN 1998S, or perhaps higher ionization at lower densities. 

Second, PTFlliqb shows an interesting and significant 
difference in its late-time Ha profile evolntion. Both objects 
exhibit qnalitatively similar asymmetric and multipeaked 
Ha lines at late times. However, PTFlliqb shifts from hav¬ 
ing a prominent bine peak at days 120-200 to a very strong 
red peak after day 500. SN 1998S does not do this, and stays 
with a bine peak continually. This is significant because the 
persistent blue peak of SN 1998S at late times was attributed 
to dust formation blocking the receding parts of the system 
llLeonard et al.ll2000l : IPozzo et af]|2004l l. The nearly identi¬ 
cal temporal evolution of PTFlliqb, but with the difference 
of ending with a stronger redshifted peak, implies that this 
asymmetry may be caused by nonaxisymmetric structure 
and viewing orientation, rather than preferential extinction 
of the far side by dust. The asymmetry in the CSM is dis¬ 
cussed more below (see §6). 

Aside from these two points, the spectra of PTFlliqb 
and SN 1998S are basically the same. How does this fit with 
their light cnrves? SN 1998S was roughly a factor of 4 more 
luminous at peak, consistent with stronger CSM interac¬ 
tion, whereas the two had very similar luminosities from 
about day 80 onward. In the next section, we show that we 
can approximate the light curve of PTFlliqb by adding the 
luminosity from wealjj CSM interaction to the light curve 
of a normal SN H-P. This veils the sharp drop in a nor¬ 
mal plateau, making it a more gradual “shoulder” in the 
light curve. If we crank up the early-time CSM interaction 
in this scenario, the peak Inminosity will be higher, and the 
shoulder in the light curve owing to the underlying SN II-P 
plateau will be harder to recognise, resembling SN 1998S. 
Thns, the inference that PTFlliqb was like SN 1998S bnt 
with weaker CSM interaction seems qnalitatively consistent 
with both the spectra and light cnrves. 


5 A MODEL LIGHT CURVE 

In Figure [S] we demonstrate that one can approximately 
match the light curve of PTFlliqb using a SN II-P light 
curve that has extra luminosity from CSM interaction added 
to it. The SN II-P light curve we adopt is that of SN 1999em 


® Here we intend “weak” to mean that CSM interaction is weak 
enough that it does not dominate the visible wavelength luminos¬ 
ity of the SN ejecta photosphere in these early phases. 


llLeonard et al.l[2002l L but of course this is just a convenient 
reference. (We could have used a different SN H-P with a 
somewhat different light curve, and then adjnsted the CSM 
parameters accordingly.) To this we add luminosity from 
CSM interaction (Lcds; the luminosity of the cold dense 
shell), calculated from a simple model with fast SN ejecta 
overtaking slower dense CSM, and t he loss in kine tic energy 
is converted to radiation (see, e.g., ISmi thl l2013all^ . Mod¬ 
els of this type are not unique and have degeneracy in ad¬ 
justable parameters like explosion energy, CSM mass, CSM 
radial distribution, and geometry. The model shown here 
is meant as a plansibility argument that PTFlliqb can be 
explained with a normal core-collapse SN that has weaker 
CSM interaction than SN 1998S, and provides only a very 
rongh estimate of the CSM conditions. 

Observations of the CDS velocity at late times when the 
emission is optically thin can help restrict some of the degen¬ 
eracy in the model. Velocities of the blue and red peaks ob¬ 
served in the late-time Ha profiles are plotted in the bottom 
panel of Figure [S] Normally, decelerating the fast 10®^ erg 
SN ejecta to the observed coasting speed of around 3000 km 
s“^ wonld require a high CSM mass, and also that roughly 
half the explosion kinetic energy was radiated away — this 
would produce a value of Lcds that is much higher than 
observed in PTFlliqb. To make a CSM-interaction model 
consistent with the observed velocities without overproduc¬ 
ing the luminosity, we must adjust the emergent luminosity 
by an artificial efficiency factor of 15%. One option to 
accomplish this physically is to simply lower the CSM den¬ 
sity and the SN ejecta mass (and SN energy) to 15%. The 
lower SN ejecta mass and energy wonld, however, then seem 
inconsistent with the underlying plateau light curve and ra¬ 
dioactive decay tail that indicate a normal SN II-P. A second 
and simpler option is to interpret the 15% efficiency as a ge¬ 
ometric effect — i.e., if the dense CSM only intercepts 15% 
of the solid angle seen by the SN ejecta, one can still have 
a normal 10®^ erg explosion and high-density CSM, but the 
resulting Lcds is much lower than it wonld be in a spheri¬ 
cally symmetric model with the same parameters. The 15% 
geometric covering factor could, in principle, arise from CSM 
that is in a disc or torus, a nonaxisymmetric shell, or clnmps 
that intercept an equivalent fraction of the solid angle. 

Agreement between our model and observations does 
not provide conclusive evidence of asymmetric CSM, bnt 
there are a number of reasons why nonspherically sym¬ 
metric CSM is plausible. First, there is clear observational 
precedent: the similar explosion SN 1998S was inferred 
to have a highly aspherical, perhap s disc-like or toroida l 
CSM based on spectropolarimetry llLeonard et al.l I2OOC1I L 
SN 1998S was also inferred to have a disc-like CSM based 
on the mnltipeaked asymmetric Ha line profiles seen at late 
times llLeonard et al.ll2000l : lGerardv et all2000l : [Fassia et al] 
I2OOII: iPozzo et al.ll2004|j . which are very similar to those 
of PTFlliqb. A flattened CSM geometry was also in¬ 
ferred for the more recent event SN 2009ip based on en- 
ergetic arguments, spectropolarim e try, and spectral clues 
JSmlHi^^laueriimi^fc Prieto|^^ et al.l l2014l : 


I Levesque et al. 20141 : Ofek et all 2013d l. Second, our spec¬ 
tra do not exhibit the narrow P-Cygni absorption that one 
might expect if the required preshock CSM density is along 
our line of sight to the SN photosphere. Last, evolved mas¬ 
sive stars commonly show highly aspherical CSM when their 
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Figure 8. A simple model for the light curve of PTFlliqb. Top panel: Luminosities based on the observed r and i-band light curves (red 
and black data, respectively) of PTFlliqb, compared to a simple model. The thick blue curve that approxi mates the light-curve shape 
(labeled LQng+1999em) is a combination of a normal SN II-P light curve (SN 1999em, shown in grey; from [Leonard et al.ll200^ and a 
simple CSM interaction model (orange curve labeled Lcds)- This model results from the radial density distribution shown in the inset 
at upper right (solid black), but where we have assumed an efficiency of 15% to account for aspherical geometry (see text for details). 
The dashed purple curve in this panel is the equivalent density distribution with the same mass if the mod el were spherically symmetric. 
The dashed green curve is the radial density distribution from the model for SN 1998S bv IChug^ ll200lf) . In the luminosity plot (top), 
the dotted blue and dotted red/orange curves show what the total luminosity would look like with a consta nt wind densi t y par ameter w, 
corresponding to the dotted black line in the radial de nsity plot. (The f?-b and light curve of SN 1998S from iFassia et ahl J 2 OO 0 II is shown 
for comparison, as are X-ray data for PTFlliqb from lOfek et al.1 ll2013bl ~)l. Bottom panel: The evolution of the cold dense shell (CDS) 
radius (black) and velocity (rust coloured) with time in the CSM interaction model from the top panel. The solid curve labeled VcDS is 
the velocity for the solid black density distribution in the upper-right inset, and the dotted orange curve corresponds to constant w for 
comparison. The blue and red data are velocities of the blue and red peaks of the asymmetric Ha line (see Figure O. 


nebulae are spatially resolved (this is discussed more in the 
next section). 

A substantial fraction of Lcds in our model should 
escape as X-rays rather than being converted to visible- 
wavelength continuum or Ha emission. Although the ob¬ 
served X-ray luminosity at early times shows that only a 
small fraction of t he X-ray luminosity escapes (Figure|8]and 
lOfek et al.ll2013l3 ). there have not yet been published con¬ 
straints for late-time X-ray emission from PTFlliqb. This is 
expected theoretically fo r most SNe Iln (IChevalier fc Irv^ 
12012 ! : I^irski et ^l2012l ). From our analysis of the late-time 


Ha (see below) we infer that the wind density parameter in 
the outer wind may actually be similar to that in the very 
inner wind, arguing that the drop in density at cm 

used to model the light curve merely reflects a decreasing 
optical depth of the outer wind and larger escape fraction 
of X-rays. As the material expands to larger radii and the 
optical depth drops, we would indeed expect a larger frac¬ 
tion of the total postshock luminosity to escape as X-rays 
(moreso if the geometry is not spherically symmetric). This 
would occur after the plateau. 

The model that matches the observed light curve best 
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in Figure [8] is qualitati vely very simi lar to the model for 
SN 1998S presented by[Chuga| ll200lll . with a higher den¬ 
sity inner wind transitioning to a lower density outer wind 
(again, this drop may merely reflect a larger escape fraction 
of X-rays at lower optical depth or a different geometric cov¬ 
ering factor; our simple model cannot constrain this further). 
The radial density profile we adopt is shown by the black 
line in the upper-right inset of Figure[Hl Note, however, that 
this is the wind density that occupies only 15% of the solid 
angle, following our assumptions about geometry discussed 
above. Therefore, this same inset to Figure [5] shows a dashed 
magenta curve, which is what the model wind density would 
be if the same mass were averaged over dvr steradians; this 
is more appropriate when comparing to the spherica lly sym¬ 
metric wind density in the model for SN 1998S by IChugaa 
(l200lll . shown by the green dashed curve in the same in¬ 
set. So, comparing these one can see that the CSM around 
PTFlliqb is qualitatively very similar to that of SN 1998S, 
but with less mass and with a division between the inner 
and outer wind at ~ cm instead of at ~ 10^® cm. 

Our adopted CSM has an effective (i.e., spherically av¬ 
eraged) wind density parameter {w = M/vw) of ui = 10^® g 
cm“^ for the inner CSM and 2.5 x 10^"^ g cm“^ for the outer 
CSM. With M = wVw, these densities translate to pre-SN 
mass-loss rates of roughly 1.5 x 10“"^ M© yr~^ (inner) and 
4 X 10~® Mq yr“^ (outer) for a wind speed of 100 km s“^. 
The total mass in the inner shell is about 10“® Mq lost in 
the ~ 8 yr before core collapse, and about 0.04 M© lost 
in the ~ 1000 yr before core collapse (out to ~ 3 x 10^^ 
cm). These are only order-of-magnitude estimates, owing to 
possible variation in wind speed, X-ray escape fraction, and 
geometry. The CSM mass is small compared to the several 
Mq in the massive expanding SN ejecta, but the CSM can 
be dense enough to resist the expansion of the SN ejecta if it 
is in a thin disc that intercepts a small fraction of the solid 
angle of the explosion. These pre-SN m ass-loss parameters 
are less extreme than for other SNe I ln jKiewe et al.ll2012 : 
ISvirski et al.l l2012l : lOfek etld] l2013bl j, and are within the 
range of mass-loss rates of several observed classes of mod¬ 
erately massive evolved supergiants (RSGs with dense winds 
like OH/IR stars, YSGs , B[e1 s uperg iants, and relatively 
low-luminosity LBVs; see lSmithI (120141 ) for a general review 
of mass-loss rates in evolved massive stars). The GSM for 
PTFlliqb does not require the most extreme levels of LBV- 
like eruptive mass loss, but it does imply some significant 
episodic ejection or wind modulation in the decade before 
core collapse. Both PTFlliqb and SN 1998S had extended 
GSM consistent with RSG winds, alt hough on the high 
end o f known mass-loss rates for RSGs dMauerhan &: SmithI 

l2012|j . 

In this model, the total energy radiated by CSM inter¬ 
action in the first 100 days is at least 10^® ergs, although the 
early rise may be caused partly to shock-breako ut luminos¬ 
ity re processed by the dense inner wind (see, e.g.. lOfek et al.l 
l2010li . This is the extra luminosity that is needed in addi¬ 
tion to a normal SN II-P light curve. This extra energy can 
be reduced somewhat if the underlying plateau SN is more 
luminous (although it can be at most a factor of ~ 1.6 more 
luminous to not exceed the observed day 80-100 luminos¬ 
ity). The amount of extra radiated energy in only the first 
20 days indicated by optical photometry is about 7 x 10“^® 
erg, but this ignores a bolometric correction and so the true 
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Figure 9. Late-time Ha luminosity of PTFlliqb, as com¬ 
pared to other SNe with strong late-time Ha from CSM in¬ 
teraction. We include SN 1998S and SN 1980K from Figure 4 
of iMauerhan fc SmithI |2012l ). with SN 1980K Ha data from 
lMilisavh evic*er'arril2 ^2h. Ha luminosities for SN 1993J are com- 
piled b\ oChandr^et^ni(20n9h . 

radiated energy in this time is probably a factor of several 
higher, given the high >25,000 K temperature at the ear¬ 
liest times. In any case, this energy is several times higher 
than the typical energy of a UV flash associated with shock 
breakout in a RSG without a dense inner wind. 


6 LATE-TIME Ha AND ASYMMETRY 

After about 100-120 days, corresponding to the drop in 
continuum luminosity of the underlying plateau SN light 
curve, the spectrum of PTFlliqb took on a different char¬ 
acter, dominated by a nebular ejecta spectrum plus a very 
strong Ha line that we attribute largely to CSM interaction. 
As noted above, the multipeaked Ha li ne profile is qualita¬ 
tively very s imilar to that of SN 1998S ([Leonard et ai.l[2000l : 
iFassia et al.ll200ll : IPozzo et al.ll2004l l. which was interpreted 
as arising from CSM interaction in a flattened disc. We note 
that the late-time Ha from PTFlliqb is also reminiscent 
of the late Ha emission from the nearby SN Ilb explosion 
SN 1993J llMatheson et al.ll2000al f9l. where the Ha was also 
attributed to CSM interaction in an extended disc. 

Figure m shows the temporal evolution of the integrated 
Ha line luminosity in PTFlliqb and in a few other SNe with 
strong late-time HaQ For the first ~ 500 days, PTFlliqb 
decays at roughly the ®®Co rate, and is quite similar to 
SN 1998S. It is during this time frame that PTFlliqb and 
SN 1998S both show a strong blueshifted peak in Ha, along 

^ Note that at late times, some of the spectra are contaminated 
by very narrow nebular Ho and [N ll] emission from adjacent 
H II regions, where variable amounts are present in the spectra 
owing to seeing differences. We subtracted the contribution of 
these narrow H II lines from the integrated line flux. 


© 2012 RAS, MNRAS 000, [T]-?? 











































16 Smith et ai 





-100 -50 0 50 100 

Velocity (km/sec) 


Figure 10. These panels relate to the discussion of asymmetry in §6 of the text. The upper-left plot compares the late-time asymmetric 
Hcf profiles of SN 1998S and PTFlliqb. The lower-left panel is a sketch of CSM interaction for a SN running into a flattened disc-like 
CSM that has a lower density on the left side and a higher density on the right. A high-inclination observer at position A would see 
stronger redshifted peaks, as in PTFlliqb, whereas an observer seeing the same event from position B would see a blueshifted peak as in 
SN 1998S. An observer at low inclination (position C) would see a narrower and more symm etric line prof i le. Th e right panels correspond 
to observations of the torus around the mass-transferring eclipsing binary RY Scuti, from [Smith et ^ (120021 '). as seen in [N ll] A6583 
emission. The botton panel shows the spatially resolved position-velocity diagram {HST/STIS spectra) with the slit aperture running 
through the major axis of the torus. This reveals an azimuthally asymmetric density distribution around the torus, with higher densities 
on the receding side. The top-right panel shows the integrated line profile from this circumstellar nebula, seen in a ground-based echelle 
spectrum. The radius of the torus is about 1000 AU, similar to the CSM encountered by PTFlliqb at late times, and the deduced 
CSM mass is similar as well. It seems likely that if a SN were to expand into a torus like that around RY Scuti, it would produce CSM 
interaction signatures much like those observed in PTFlliqb at late times. 


with strong nebular lines like the broad Ca ll IR triplet. Dur¬ 
ing this time, the Ha luminosity may be powered largely by 
radioactive decay. As such, the velocity asymmetry in the 
Ha lines (especially the blue bump), could indicate a highly 
nonspherical distribution of ®®Ni. 

After about 500 days, however, the Ha from PTFlliqb 
decays more slowly than in SN 1998S, SN 1993J, and 
SN 1980K, and no longer traces the ®®Co radioactive-decay 
rate. During these later times, Ha is therefore most likely 
powered by CSM interaction. It is during these same later 
epochs after day 500 when PTFlliqb’s Ha profile switches 
from a blueshifted peak to a stronger redshifted peak. It 
will be interesting to see how PTFlliqb continues to evolve. 
In SN 1998S, the drop in Ha luminosity reached a floor 
at around day 1000, after which it followed a trajectory 
consistent with CSM inter action in a very dense RSG wind 
dMauerhan fc Srn ithll201 2l). a nd SN 1980K and SN 1993J be- 
haved similarly jFesm et al.lll99^ : iMilisavlievic et aLll2012l : 
iMatheson et al.ll2000bl) . The stronger day 500-1100 Ha lu¬ 
minosity in PTFlliqb could indicate an earlier onset of this 
behaviour in a wind that is an order of magnitude more 


dense (which would be extreme, given that the wind pa¬ 
rameters for SN 1998S resembled those of VY CMa, the 
strongest known RSG wind), or perhaps it is merely over¬ 
taking a shell or terminal shock in the RSG wind that tem¬ 
porarily enhances the strength of Ha. Continued observa¬ 
tions of PTFlliqb are encouraged. In any case, the fact that 
PTFlliqb had weaker CSM interaction at early times and 
stronger interaction at late times compared to SN 1998S sug¬ 
gests that the pre-SN mass loss in SNe Hn does not necessar¬ 
ily follow a common recipe, and that detailed observations 
of the shock progressing through the distant CSM may help 
unravel their highly varied mass-loss histories. Most SNe are 
not observed sufficiently late to see this behaviour, and/or 
they are too distant. If pre-SN mass loss is punctuated by 
stochastic mass ejections over a range of timescales, this 
would be quite relevant for diagnosing the driving mecha¬ 
nism in the final phases before core collapse. 

Similarly, if the stochastic mass loss is highly asym¬ 
metric (especially azimuthally asymmetric), this may also 
be an important clue to its underlying physical mechanism. 
When the asymmetric blueshifted Ha line profiles were seen 
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at late times in SN 1998S, it was inferred that dust ob¬ 
scuration may have played an important role in making an 
otherwise more symmetric double-peaked line appear he^- 
ily blueshifted ll/eonard et alJ [ioM: iGerardv et ^ I 2 OOOI : 
IPozzo et al.ll2004l ~l. IPozzo et alJ ( 2004l i mentioned the pos¬ 
sibility that dust formed in the postshock region (the CDS) 
may have been contributing to the asymmetry, as was seen 
later i n the unambiguous case of SN 2006jc dSmith et alJ 
l2008bh . However, as noted earlier, we cannot rely upon dust 
or ejecta obscuration to account for the strong asymmetry 
seen in PTFlliqb, because this time it is the blueshifted ma¬ 
terial on the near side that is missing. The fact that emission 
from the far side of the CSM interaction is brighter provides 
a strict requirement that the density distribution in the CSM 
is intrinsically asymmetric. Interestingly, the velocity of the 
brighter red peak is slower (roughly -1-1100 km s“^) as com¬ 
pared to that of the blue peak (roughly —2000 to —3000 km 


denser CSM is an intuitively plausible explanation for the 
brightening of the red side of the line. If SNe Iln typically 
have such asymmetry, then a range of viewing angles could 
easily explain why some have redshifted peaks, some have 
blueshifted peaks, and others have narrower and centrally 
symmetric Ha profiles at late times, as depicted in the left 
panels of Figure [TOI 


What sort of process could produce a flattened CSM 
geometry that is significantly more dense on one side than 
the other? Axisymmetric structures in the CSM (discs and 
bipolar nebulae) may conceivably arise from rapid rotation 
in single stars, but high degrees of azimuthal asymmetry are 
difficult for single stars to achieve. A way to have one-sided 
CSM while still satisfying the equatorial mass distribution 
suggested by spectropolarimetry of SN 1998S is to have mass 
loss in a binary with some nonzero eccentricity, or with un¬ 
steady mass loss. As an illustrative example, the right panels 
of Figure [To] show observed emission from the nebula around 
the Galactic massive star RY Scuti. RY Scuti is a rare exam¬ 
ple of a massive eclipsing binary system caught in the brief 
phase of mass transfer, where one star is being stripped of 
its H envelope on its way to the WR phase, and it is so 
far the only such sys tem known with a spatially resolved 
toroidal CSM nebula (jCrundstrom et ^l2007l : ISmith et al.l 
l2002l . [2C)llbl f . Figure [To] shows the spatially resolved velocity 
structure of the ring and the integrated emission-line profile 
of [N ll] A6583, from ISmith et af]l2002l '). The integrated line 
profile from the nebula is multipeaked and very asymmetric, 
with a brighter red peak — very much like PTFlliqb at late 
times, but with slower (preshock CSM) expansion speeds of 
only ±40 km s“^. When emission from RY Scuti’s torus 
is spatially and spectrally resolved (see the position-velocity 
diagram in the bottom right in Figure fTOll . it is clear that the 
density distribution in the expanding torus is azimuthally 
asymmetric, with several clumps ar ound the torus an d gen¬ 
erally higher density on the far side ijSmith et al.ll200^ f. This 
azimuthal asymmetry arises despite the fact that the cen¬ 
tral eclipsing binary system appears to be tidally lo cked 
and has a circularised orbit (jCrundstrom et aP 120071 1 . In 
this case, the azimuthal asymmetry may arise from episodic 
mass loss. Proper motions of the multiple shells around 
RY Scuti r eveal two separate major ejections in the last 
~ 250 yr ([Smith et al.l boilbll . Such sequential mass-loss 
episodes in a binary might provide the complex density 



Figure 11. Early-time spectra of PTFlliqb (day 2 after discov¬ 
ery, roughly day 5 af ter explosion), SN 20 13cu (day 3 after dis¬ 
covery and explosion: iGal-Yam et al.l l201^h_ajid_SN 1998S (day 4 
after discovery, day 5 after explosion: [Leonard et aLll2000h show¬ 
ing the WR features. The days noted next to each SN name are 
approximate times after the inferred explosion time. Various lines 
of highly ionised species such as He I, He II, N III, and C III are 
seen in the spectra, plus broad Lorentzian wings of Balmer lines. 
These WR-like features appear despite the likelihood that both 
the PTFlliqb and SN 1998S progenitors were RSGs, not WR 
stars, because the dense cool wind becomes highly ionised by 
shock breakout or CSM interaction from the fastest ejecta. 


structure around PTFlliqb; interestingly , the total mass in - 
ferred for RY Scuti’s nebula of 0.001 Mq ([Smith et al.ll200^ 
within 1000 AU is in the right ballpark. We infer that if a 
Type II-P explosion was surrounded by RY Scuti’s nebula, 
the resulting CSM interaction could produce the late-time 
Ha emission seen in PTFlliqb. The analogy can only be 
pushed so far, however, because RY Scuti is a binary of an 
O star and a B supergiant, so it is unlikely to produce a 
plateau light curve if it exploded tomorrow. Nevertheless, 
it provides a nice illustration of the asymmetry that can 
arise in the CSM around interacting binaries. Many other 
possibilities exist to create azimuthally asymmetric struc¬ 
ture, such as wind collisions or other interactions in binary 
systems with nonzero eccentricity. 


7 WOLF-RAYET FEATURES IN THE EARLY 
SPECTRUM 

PTFlliqb now joins a few other SNe (Types Hn and Hb), 
which have been caught early and show the WR features 
associated with He ll A4686 and strong Lorentzian wings 
of Balmer lines in their very early-time spe ctra. The best- 
studied cases of this so far include SN 1998S ( Leonard_eLalJ 
200^ Chugaj] I 2 OOII ) and SN 2013cu JCal-Yam et al.l 20141 : 
Groh[ 2 OI 4 I I. Figure [TT] shows these features in detail as 


they appear in the early spectra of PTFlliqb, SN 1998S, 
and SN 2013cu. The spectrum of SN 2013cu from 
[Cal-Yam el^jjj __([20]T 1_ was obtained from the WISeREP 
database ( Yaron fc Cal-Yamll2ni2l l. 
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Recently. I Gal-Yam et al.l il2014l ~) reported the detection 
of these same features in the first few days after explo¬ 
sion in SN 2013cu, which was a stripped-envelope Type lib 
event. Those authors interpreted the spectral features as re- 
semblinK a WN6h sp ectral type. Based on this similarity, 
iGal-Yam et al.l (l20l4) proposed that the progenitor was a 
WR-like star. This is an important claim, since no WR pro¬ 
genitors have been detected yet for any stripped-envelope 
SNe, but relatively cool YSG progenitors have been di¬ 
rectly detected in other SNe Ilb, the thre e clear cases being 
SN 1 993J, SN 2011dh, an d SN 2013df llMaund &: Smar^ 
I 2 OO 9 I : IVan Dvk et al.ll2013l . l2014h . (Interestingly, SN 1993J 
behaves very similarly to PTFlliqb and SN 1998S at late 
times.) 

However, in this paper we have shown that these same 
WR features are seen in PTFlliqb, and they had also been 
seen previously in SN 1998S|f] This chall enges the WR- 
like p rogenitor interpretation forwarded by iGal-Yam et al.l 
I 2014h. since both P TFll iqb (this work) and SN 1998S 
i Mauerhan &: Smithll201^ : IShivvers et al]l20l4 ) have GSM 
consistent with RSG winds. In the particular case of 
PTFlliqb, we see a plateau in the light curve, indicating 
that the underlying SN event was a Type II-P explosion — 
this requires that the progenitor was an extended cool su¬ 
pergiant at the moment it exploded, and not a compact WR 
star. From this comparison alone we cannot rule out the pos¬ 
sibility that SN 2013cu had a WR-like progenitor, but clearly 
a WR progenitor is not necessary to yield WR features in a 
SN. A YSG progenitor, as has been detected directly in other 
SNe lib, may provide a suitable alternative for SN 2013cu. 
Indeed, compari ng the spect rum of SN 2013cu to radiative- 
transfer models. ICrohl ||2014| ') concludes that the SN 2013cu 
progenitor could not have been a WR star, and instead 
must have been a yellow hypergiant (with more extreme 
mass loss than a normal YSG). The enhanceme nt of N lines 
in the spectrum that led iGal-Yam et al.l ll2014l ~) to connect 
SN 2013cu with a WN star (rather than a WC) star is not 
exclusive to WR stars, since many cooler evolved massive 
stars across the upper HR diagram show N-enriched atmo¬ 
spheres (in any case, the strong N lines may b e more of a 
temperature effect than abundance; iGro Eioli). Neverthe¬ 
less, early time data such as the spectra obtained for 2013cu 
and PTFlliqb pr ovide powerful prob es of the immediate 
pre-SN mass loss (iGal-Yam et al.|[2014ll . 

We conclude that the WR-like spectrum seen at early 
times has little to do with the spectral type of the progen¬ 
itor before explosion. Instead, it most likely reflects a wind 
density that is high enough to be optically thick at radii of 
~ 10 AU after it becomes ionisedlj Interestingly, the im- 


® We note one caveat, related to sensitivit y or time resolution. An 
emission feature at ~7100 A identified by iGal-Yam et al.l (l2014l i 
as N IV was strong in the spectrum of SN 2013cu taken 15.5 hr af¬ 
ter explosion, but the feature faded beyond detectability by day 3 
(see Fig ureLUJ. This feature was also seen in SN 1998S, although 
[Leonard et al.l ll200Clh identified it as a C II blend. Unfortunately, 
our earliest spectrum of PTFlliqb is too noisy at these waye- 
lengths to see if the feature is present at the same level as in 
SN 1998S at a similar time. 

® The estimate of ~ 10 AU in the case of PTFlliqb comes from 
the required photospheric radius, given the observed luminosity 
and temperature we quote here. A radius of 10 AU is also consis- 


plied wind den sity requires a mass-loss rate (0.0 01 Mq yr“^ 
for SN 2013cu; iGrohll^Oll : iGal-Yam et al.ir2014l l that is or¬ 
ders of magnitude higher than observed winds of WR stars, 
and well above the line-driving limit (see lSmi ^ l2014ll : WR 
stars do not have such high densities at 10 AU because their 
winds are fast. The required density is more easily achieved 
in the slow dense winds of cool supergiants (with lower wind 
speeds, a lower mass-loss rate is required for the same w). Of 
course, cool supergiants normally lack the ionising photons 
needed to create the characteristic WR spectrum. A WNH- 
like spectrum will be produced, however, when any dense 
cool wind is bathed in a huge flux of hard ionising photons 
that are capable of fully ionising the wind. The WN features 
can be strong even when H is abundant in the atmosphere; 
as a consequence, there are many WNH stars that are not 
at the same evolutionary stage as He -burning WR stars (see 
ISmith fc Contiir2008l : lGrowtherll2007l ). The requisite blast of 
hard photons might come f rom either the UV fl ash asso¬ 
ciated with shock breakout (lGaTYam_et_aL 20l4l. or from 
GSM interaction ([Leonard et al.l 120001 : IChusail 1200 il l gener¬ 
ated when the fastest outer SN ejecta begin to crash into 
the slow inner GSM. Since this occurs inside the continuum 
photosphere, it is difficult to know which process dominates, 
but nevertheless, either process requires a similarly optically 
thick inner wind. 

Lastly, we note that independently, Shivvers et al. 
(2014) draw a similar conc lusion regarding the interpreta- 
tion of the WR features bv iGal-Yam et ^ (|2014| L Shivvers 
et al. present a high-resolution echelle spectrum of SN 1998S 
obtained at very early times when the WR emission features 
were seen in that object. They find a slow wind consistent 
with a RSG progenitor, and therefore favour a similar inter¬ 
pretation that a WR-like spectrum at early times does not 
necessarily indicate a WR progenitor. 


8 OVERVIEW OF PTFlliqb: BRIDGING TYPE 

Iln EVENTS WITH NORMAL SNe 

In principle, any type of SN explosion could appear as a 
Type Hn, since the Hn spectroscopic signature is not the re¬ 
sult of the explosion physics, but is caused instead by H-rich 
(or H-poor/He-rich in the case of SNe Ibn) GSM interaction 
that happens afterward. When GSM interaction is strong, 
the GSM optical depth is high. Gonsequently, in cases where 
GSM interaction substantially enhances the peak visible con¬ 
tinuum luminosity of a SN Hn, the photosphere can actu¬ 
ally be ahead of the forward shock in the GSM, and hence, it 
can effectively mask the identity of the underlying SN ejecta 
spectrum (see, e.g.. [Smith et al.1[2008a[ for more details). 

Because of its relatively weak GSM interaction at early 
times, PTFlliqb affords an interesting opportunity to more 
clearly see the underlying SN H-P photosphere. The signa¬ 
ture of a SN H-P is seen from the shoulder in the light curve, 
which can be explained as a composite of a normal plateau 
and GSM interaction, and the visible-wavelength spectrum 
toward the end of that plateau matches a SN ILL, when the 


tent wi th the radius of the C DS on day 2 in our CSM interaction 
model. ICal-Yam et al.l |[2014l ) quote a similar approximate radius 
of 10^^ cm. 
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broad-lined SN ejecta dominate the observed light. More¬ 
over, it seems clear from the composite SN II-P plus CSM 
interaction light-curve model that all one needs to do is 
“crank up” the CSM interaction by a factor of a few (in 
CSM density and luminosity) to make this object look like 
the more traditional SN Iln explosion SN 1998S. 

When we first described the spectral evolution of 
PTFlliqb in §3.2, we highlighted three main phases; 

(1) The early phase (days 0-20) dominated by a 
smooth blue continuum, narrow emission lines with broad 
Lorentzian wings, and the WR-like spectral features. During 
this time, the spectrum was that of an opaque wind ionised 
by either CSM interaction or a UV flash from shock break¬ 
out (or both). 

(2) The plateau phase (days 20-120) when the spec¬ 
trum resembled the broad-lined photospheric spectrum of a 
normal SN II-P or II-L, with broad emission and absorption 
profiles. In this phase, the emitted light is dominated by the 
fast SN ejecta, as in a normal SN. 

(3) Late-time CSM interaction with strong asymmetric 
Ha, which can be subdivided into the early nebular phase 
(days 120-200) when Ha had a strong blueshifted peak and 
was probably powered largely by radioactive decay, and later 
times (after day 200) when Ha must have been powered by 
CSM interaction and had a stronger redshifted peak. 

In Figure 1121 we show a sketch that illustrates a pro¬ 
posed sequence of events to qualitatively explain this ob¬ 
served evolution of PTFlliqb. The four panels in this figure 
are described below. 

Figure I12h : The progenitor is a RSG or yellow super¬ 
giant/hypergiant surrounded by a dense wind and an even 
denser disc-like distribution of CSM. The densest parts of 
the disc may reside at a radius of ~ 10 AU. The disc is likely 
the remnant of some previous binary interaction events. 

Figure 112b : For the first 20 days immediately after 
explosion, the SN ejecta crash into the inner dense wind of 
the progenitor at a radius of ~ 10 AU. However, the wind 
is opaque once it is ionised, and a photoionised precursor 
(due to either CSM interaction X-rays that are reabsorbed, 
or to a UV flash from shock breakout) creates an electron¬ 
scattering photosphere outside the fast SN ejecta. Narrow 
lines are emitted by the preshock wind. The collision in the 
equator slows the expansion there, but the fast ejecta expand 
freely at other latitudes. 

Figure ll2b : At days 20-120, the radius of the SN ejecta 
photosphere has now grown to ~ 100 AU at latitudes away 
from the equatorial plane, as in a normal SN II-P. Although 
the photosphere is constantly receding in mass coordinates 
through the expanding ejecta, it is larger than the radius 
where the most intense CSM interaction is occurring in the 
equator (10-20 AU). Thus, the equatorial CSM interaction 
region is enveloped hy and hidden inside the normal broad- 
lined SN photosphere as seen from most viewing angles. This 
is the most likely explanation for the disappearing or weak¬ 
ening narrow lines in some SNe Hn. If CSM interaction still 
contributes to the total plateau luminosity, the photons must 
diffuse out through the fast ejecta. 

Figure 12. Illustration of a proposed sequence of events to ac- Figure [T2ji: At late times after the plateau drop (days 

count for PTFlliqb’s evolution. The four panels la—d) are dis- .i .• i ^ i i j. n 

, . , ' ' 100-120), the continuum photosphere recedes to smaller 

cussed m the text, 58. i • 

radii and once again exposes the CSM interaction occur¬ 
ring in the disc. As the photosphere recedes completely and 
the ejecta become optically thin, the emergent spectrum is 



© 2012 RAS, MNRAS 000, [T]-?? 
































20 Smith et al. 


dominated by optically thin ejecta powered by radioactive 
decay plus CSM interaction in a disc or torus. The postshock 
speed of the swept-up CSM is a few 1000 km s“^, indicated 
by the width of the multipeaked asymmetric Ha line. 

This sequence of events gives a satisfactory explana¬ 
tion for the spectroscopic evolution seen in PTFlliqb. By 
changing the density and radial extent of the CSM disc (and 
perhaps also by changing the viewing angle of the observer) 
one can adapt this scenario to other SNe; with denser or 
more extended CSM, we get a SN Iln like SN 1998S, and 
with less dense or less extended CSM, we may get normal 
SNe II-L or II-P events. Thus, this scenario for PTFlliqb 
provides a direct bridge between SNe Iln and SNe II-P/L. 

This, of course, does not mean that all SNe Iln are 
actually Type ITP explosions underneath their CSM inter¬ 
action, but it does confirm that some of them are under¬ 
lying SNe II-P that arise from RSGs or YSGs, although 
masked at some phases by GSM interaction. Most of the 
interest in SN Iln progenitors has thus far concentrated on 
LBVs, partly because in some cases the GSM mass is so high, 
and also because the very bright LBV-like progenitors may 


be easier to detect in pre-explosion data (iGal-Yam et al 
200^; [GahTam&Leonar^ [200^ SmithI [20141: ISmith et al 


20071. l2008al . l2010bl. l2011cll . ISmith et all (^2009a^ suggested 


that extreme RSGs like VY CMa have CSM environ- 
ments that can produce lower-luminosity SNe Iln, and 
iMauerhan fc SmithI ll2012ll showed that the late-time CSM 
interaction in SN 1998S was consistent with such a wind 
operating for ~ 1000 yr before core collapse. More recently, 
IShivvers et al.l (l2014h have also found evidence for a RSG 
progenitor for SN 1998S based on very early high-resolution 
spectra, as noted above. There have been other suggestions 
of RSG progenitors of some lower-luminosit y SNe Iln as well 
llStritziriger et af]l2012l : ISmith et al.ll2009^ . 

From the point of view of understanding the connection 
between SNe and their progenitors, an important question 
then arises: For how long before explosion was PTFlliqb’s 
progenitor a cool supergiant? Extreme instabilities in the 
years leading up to core collapse are required to eject many 
Mq in some cases that produce luminous SNe Iln, but there 
may also be less extreme instabilities that could drive a 
wind, cause violent pulsations, or simply inflate the star’s 
envel o pe in the final year s of a s tar’s life dSmith _|^ Arngt3 
20141: lo uataert fc Shiod^ l2012l: IShiode fc Quataerd 12014 


Qfek et al. 2014ah . Yoon fc Gantiell3 ~i 2010ll have also sug¬ 
gested that some luminous RSGs may experience an en¬ 
hanced superwind phase for ~ 1000 yr before core collapse 
driven by pulsational instabilities. During most of its He- 
burning lifetime, the progenitor may therefore have looked 
quite different from how it looked in the decades or cen¬ 
turies immediately preceding explosion. Whatever the phys¬ 
ical cause, it does seem clear that PTFlliqb resides between 
SN 1998S and normal SNe II-P in terms of its luminosity and 
CSM interaction intensity. 

Thus, this type of pre-SN activity, instability, and vari¬ 
ation in stellar structure may play an important role in un¬ 
derstanding the connection between the most common core¬ 
collapse events (normal SNe II-P) and other objects like 
SNe II-L and SNe Iln that appear to have some overlap. 
It could be that SNe II-L and SNe Iln form a sequence of 
increasin g progenitor mass and i ncreased pre-SN instability 
(see, e.g.. lYoon fc Cantielldl201^ L or perhaps it could be due 


to other factors, such as proximity of a com panion star in bi¬ 
nary systems that suffer these instabilities (ISmith fc Arnetj 
l2014h . Exploring the diversity in CSM density and geometry 
may be extremely important for sorting this out, since disc¬ 
like CSM around a RSG is unlikely without a companion 
star to supply angular momentum. As noted above, some 
objects classified as SNe II-L, including c lassic objects such 
as SN 1980K, SN 1979C , and others Isee lFesen et al.l[19991 : 
iMilisavlievic et ^l2012l L do not have spectra obtained at 
very early times — perhaps they would have been classified 
as SNe Iln if they had been caught sufficiently early. Detect¬ 
ing additional cases like PTFlliqb and understanding how 
common they are may help flesh out the continuum of diver¬ 
sity in pre-SN mass loss. If there are additional weak SNe Iln 
that are not easily recognised because the SN Iln signatures 
are fleeting, it would i mply that pre-SN instability is not 
limited to only 8-9% (ISmith et al.ll2011al l of core-collapse 
SNe that have been categorised as Type Iln. 
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